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Q ■ Abstract 
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f^ . We compute the third-order correction to the heavy-quark current correlation func- 

tion due to the emission and absorption of an ultrasoft gluon. Our result supplies 
a missing contribution to top-quark pair production near threshold and the deter- 
mination of the bottom quark mass from QCD sum rules. 
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In a previous paper [1] we described the calculation of the ultrasoft gluon contribution 
to the correlation function of heavy-quark vector currents, which constitutes one of the 
major missing pieces in perturbative calculations of quarkonium-like systems at the third 
order in non-relativistic perturbation theory. That paper presented the result for the 
residues of the bound-state poles of the correlation function, which relate to quarkonium 
decay constants. (The ultrasoft contribution to the S'-wave quarkonium masses has been 
known for some time [2j.) However, some of the most interesting physics related to the 
heavy-quark threshold - the determination of the bottom quark mass from sum rules 
[3] and the top-quark pair-production cross section [4j - requires the calculation of the 
full energy-dependent correlation function (see also the reviews 05] )• In this Letter we 
provide the result for the ultrasoft correction to the full correlation function. 

The effective field theory formalism and technical set-up for this calculation have 
already been given in [1] and will not be repeated in detail here. In brief, we consider 
the current correlation function 



2{d - l)N,G{E) = i / d'^xe'^'' (0|T[xVV](x) [V^V*x](0)|0) , (1) 

where E = ^/s — 2m is small, y/s is the centre-of-mass energy and m the heavy-quark 
pole mass, (d = 4 — 2e is the space-time dimension in dimensional regularization.) 
The heavy-quark current ip^a'^x is defined in non-relativistic QCD. The (leading) ultra- 
soft contribution to G{E) refers to diagrams where one gluon has ultrasoft momentum 
k r^ E r^ mv"^, while an infinite number of potential (Coulomb) gluons can be exchanged 
between the heavy quarks, promoting the heavy-quark propagators to the Green func- 
tion of the Schrodinger operator including the colour Coulomb potential. Computing 
Feynman integrals with Coulomb Green functions while simultaneously regulating all 
divergences dimensionally to be consistent with fixed-order matching calculations (de- 
fined according to the threshold expansion [7]) is the main challenge of the ultrasoft 
calculation. From the ultrasoft correction to ([T]), 5^'^G{E) (Eqns. (8), (9) of [1]), the cor- 
responding correction to the inclusive heavy-quark production cross section is computed 
as 

[i?]„, = 127re^-^-Im5"^G(E), (2) 

where R = (Jqqx/'^o with o"o = 4:7ial^/{3s) is the usual i?-ratio, and eg the heavy-quark 
electric charge in units of the positron charge. 

The calculation of the ultrasoft contribution to the current-correlation function in- 
volves ultraviolet divergences of various kinds. Divergences in box-type subdiagrams that 
do not contain any of the two current operator insertions are subtracted by counterterms 
related to the renormalization of the potentials in the effective Lagrangian. Another class 
of divergences arises from vertex subdiagrams with up to three loops and lines connect- 
ing to one of the current insertions; they are cancelled by the counterterms belonging 
to the non-relativistic heavy-quark current operators. These divergences have the same 
structure in the correlation function and the bound-state residues, and have already 
been discussed in [T]. In addition G{E) has overall divergences from propagator-type 
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Figure 1: Three- loop diagrams generating an overall divergence proportional to o?^Eje. 
Thick lines denote an unstable top-quark propagator, wavy (dashed) lines ultrasoft (po- 
tential) gluons; the square the non-relativistic current insertion. Symmetric versions of 
the last two diagrams are not displayed. 



diagrams (involving both operator insertions) with up to five loops, which are not rele- 
vant to the bound-state pole parameters. After subtraction of sub divergences the overall 
divergences are polynomial in the external "momentum" E and should not contribute 
to the imaginary part ([2]). However, there is a subtlety for top-quark pair production, 
which we now explain (see a related discussion in [B]). 

The issue for top quarks is its large decay width Fj ~ 1.4 GeV, which happens to 
be of order of the relevant non-relativistic energies E. The correct prescription for the 
calculation of the ultrasoft correction is to replace E' ^ E' + zFt in all equations [1]. 
Equivalently, we may consider E to be complex with a finite imaginary part rather than 
the -|-ie-prescription for stable quarks. Now, the ultrasoft calculation yields an overall 
divergence (arising from quadratically divergent three-momentum integrals) 

[(5"^G(i^)]overall OC ^ ■ E (3) 

from the three-loop diagrams shown in Figure [H If i? is complex, this results in a di- 
vergence a^Ft/e in the heavy-quark production cross section ([2]) that is not cancelled 
by any counterterm associated with the heavy-quark currents. A similar overall diver- 
gence has already appeared in NNLO calculations of top-quark pair production such as 
[9], where it arises from two-loop diagrams with both loops in the potential region as 
shown in the first diagram in Figure [2l The origin and cancellation of such divergences 
is best understood in unstable-particle effective field theory fTD], which provides a con- 





Figure 2: Cancellation of the overall divergence agE/e at NNLO against a partially 
"non-resonant" electroweak contribution involving a Wh top-quark self energy insertion. 
Symmetric versions of the last diagram are not displayed. 



sistent treatment of finite-width effects beyond the E —>■ E + iTt replacement (valid 
only up to NLO), and includes non- resonant contributions to physical cross sections, 
here e"*"e~ —^ W~^W~bb, which do not include the unstable top quark in the final state. 
In this framework, the overall ultraviolet divergence a^Ft/e from the first diagram in 
Figure 2 is cancelled against an infrared divergence from the second diagram. Note that 
the second diagram is not present in non-relativistic QCD or for stable quarks, since all 
loops are in the hard region, and that the electroweak self-energy insertion is computed 
in conventional perturbation theory in the full standard model. The second diagram 
is O(asa^), so the parameter dependence matches, since Fj ex a^, where a^ denotes 
the SU(2) gauge coupling. A similar cancellation between the overall divergences from 
Figure [1] with non-resonant contributions is expected at NNNLO. 

The ultrasoft contribution to the NRQCD correlation function provides only part of 
the third-order result for the heavy-quark pair production cross section near threshold, 
and is regularization- and scheme-dependent. Our conventions for dimensional regular- 
ization and the counterterms have been specified in [IJ and are chosen to conform with 
standard conventions for the calculation of MS-subtracted coefficient functions. With 
respect to the overall divergences discussed above, we note that we perform the calcula- 
tion of the correlation function ([1]) with NRQCD for stable quarks in the complex energy 
plane at E + iTf. This corresponds to keeping only the leading width- correct ion in the ef- 
fective Lagrangian for unstable quarks, (iFt/2) ip'^ipy i^i the framework [IQJ. The 1/e poles 
are then subtracted minimally (MS). The resulting scheme and scale dependence cancels 
with other NRQCD contributions to the correlation function, and matching coefficients, 
but there is a left-over dependence proportional to F^ due to the overall divergences 
discussed above, which cancels with electroweak corrections that are not yet known. 

We are now in the position to present our main result. (The technical details of the 
calculation will be explained together with those of [I] in a separate publication.) We 
express this in terms of the dimensionless complex energy variable E = E/{ma1), where 
Im. E = Tt = Ft/(ma^), and 

A = ^=, p = ln(^)+7^ + V^(l-A), (4) 

-E 



with Cp = 4/3, and ip{z) = dT{z)/dz the ■j/'-function {^^'{z) denotes its first derivative). 
Gc = (2/3)(l — 1/A — 2V), used below, is related to the Coulomb Green function at zero 
radial distance by Gc = m'^as/{4:7c) Gc- The imaginary part of the correlation function 
at complex energy involves divergent, logarithmic and finite contributions. We obtained 
the divergent and logarithmic terms in closed form and the remainder numericallylj We 
thus obtain 
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"^ Dropping the divergent 1/e pole terms in ([5]) amounts to MS renornialization. 
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The double logarithmic and 1/e pole terms are identical to those that appear in the result 
for the wave function at the origin (Eqn. (13) of PJ) under the replacements Tf —>■ 0, 
Gc ^ K = (8/9) {al/7r) \ij^{0)\y {2m^a{/{97r^)). By expanding ^"^^(E) around the 
bound-state poles A = n, we also reproduce the single logarithmic terms in [Ijo Only 
the imaginary part of 6'^'^G{E) is needed for the cross section ([2]), therefore in writing ([5j) 
we already dropped some purely real terms. The imaginary part of the non-logarithmic 
correction S'^'^{E) is tabulated in Tabled] for a set of values of the real part (rows) and 
imaginary part (columns) of ii^ in a range relevant to top-quark pair production. The 
table can be used to generate an interpolating function with an accuracy better than 0.1 
in the entire range of the tablejj 

The new ultrasoft correction has a very large effect on the ti cross section near 
threshold as illustrated in Figure [31 Here we adopt the top quark pole mass m,t = 
175 GeV, and fix a^ = 0.14, which corresponds to the Bohr radius scale fis = 32.5 GeV. 
The solid line in the upper panel of Figure [3] shows the non-logarithmic contribution 
from Itsi6'^^{E) to [R]us alone, which is seen to be around +25% in the peak region, in 
nice agreement with the estimate from the wave- function calculation \T\. Including the 
logarithmic term requires a choice of the scale /i. Since /i is not related to the scale 
of as, but designates a factorization scale that separates the ultrasoft from hard and 
potential contributions, we choose the two values fiB and rrit to represent a reasonable 
range. Since the factorization scale dependence is sizable this results in a large range of 
[R]us reflected in the two dashed curves in Figure [3] (upper panel). We add these two 
results to the leading order ti cross section in the lower panel. Our results show that 
despite the large quark mass, third-order perturbative corrections from the ultrasoft scale 



^To this end write 
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The correction to |?/'„(0)p is then given by (3X/4) (fe„ + 3a„/n). 

^For instance, using Mathematica's built-in function Interpolation with default setting. 
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Table 1: Value of the non-logarithmic ultrasoft correction Im5"*(i?) for various scaled 
energies E = E/{ma1). Columns refer to different YieE from —4.1 to 2.7, rows to five 
values oilmE = Tt = Tt/{ma1) for given real part. 



can have a significantly larger impact on top-quark pair production than anticipated. 
However, it should be kept in mind that the ultrasoft correction is not physical by itself 
as is clear from its factorization scale dependence. (The non-logarithmic term is still 
factorization-scheme dependent.) In [TT] we combined the ultrasoft correction reported 
here with the third-order potential correction |i2j| and all other calculated third-order 
terms and observed a significant cancellation between the ultrasoft and potential terms. 
A final assessment of theoretical uncertainties should therefore be attempted only once 
all third-order corrections to the ti cross section have been assembled. 
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Figure 3: Top panel: Ultrasoft correction [R]us to ti production as function of E = 
^/s — 2mt. Solid: non- logarithmic contribution only. Dashed: total contribution with 
/^ = A^B = 32.6 GeV (upper dashed) and n = mt = 175 GeV (lower dashed). Parameters: 
rrit = 175 GeV, Tf = 1.4 GeV, a^ = 0.14. Bottom panel: Ultrasoft correction added to 
the leading order (LO) cross section, i.e. [-R]lo + [R]us- The band is obtained by varying 
the scale /i between fj,B (upper line) to rrit (lower line). The solid black line refers to the 
LO cross section. 
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Table 2: Values of the non-logarithmic ultrasoft correction Im [5'^^{E + ie)] in the con- 
tinuum (real positive E). 



The heavy-quark correlation function near threshold is also a crucial input to the 
determination of uii, from large-moment bottomonium sum rules p^. The rath moment 
of the hh production cross section is defined by 
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Taking large n, typically n > 4, enhances the sensitivity to mi, and the threshold region, 
but requires a non-relativistic treatment to sum corrections of order as\fn to all orders 
in perturbation theory. On the other hand n cannot be too large, since mi,/n, the 
typical non-relativistic energy of a hh pair contributing to the nth moment must be 
larger than Aqcd to justify a perturbative computation [H]. To calculate the derivatives 
of the QCD vector current correlation function n(g^), we evaluate the moment integral 
in ([7]) with Ri,i{s) expressed as a sum over Coulomb resonances and a continuum for 
E = y/s — 2mb > 0, which should be dual to the corresponding integrated physical hh 
cross section. 

The ultrasoft correction [R]us to the continuum cross section is obtained from (^ 
and 

lm6'''G{E)= Mm ImS'^'GiE + tV) {E>0). (8) 

While we have an analytic representation of the logarithmic contributions, see (^, our 
numerical implementation of the non- logarithmic term S'^^{E) does not allow us to make 
the imaginary part of E arbitrarily small. We therefore calculate 5'^'^G{E + iV) for 
several values of F (for given E) and obtain the continuum value by extrapolating to 
r = through a polynomial fit. By variations of this procedure we estimate that the 
relative error in the calculation of the ultrasoft contribution to the continuum cross 
section (real i? > 0) is less than 1%. Numbers for the non-logarithmic term are provided 
in Table El 

The sum of the leading-order and ultrasoft contribution to the moments is given in 
Table E] together with the leading-order one alone for &-quark pole mass m^ = 4.95 GeV 
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Table 3: Moments of the hb spectral function for mf, = 4.95 GeV. (The moment integral is 
cut off at E = 20.) The renormalization/factorization scale is taken to be Hn = Irn}^! ^fn 
corresponding to ^^(/i^) = 0.228,0.250,0.267 for n = 6, 10, 14. The last row shows the 
sum of the leading-order moments and the ultrasoft contribution. The leading order 
alone is given in the first row for comparison [Hj. 



and renormalization/factorization scale /i„ = 2r7ib/y^. Given the large size of the ul- 
trasoft term in the tt cross section it may not be surprising that here we find that the 
ultrasoft correction is 30% to 80% of the leading-order term, putting the perturbative 
approach into doubt for the larger moments. (The correction from the non-logarithmic 
term alone is even larger, cf. Figure [3l upper panel.) However, as mentioned above, 
the result for the ultrasoft correction alone should be regarded with caution, and there 
is reason for assuming that the large ultrasoft correction is partially a consequence of 
MS factorization, such that the true third-order correction is smaller when all other 
corrections are added. 

To summarize, we evaluated the correction to the vector-current heavy-quark cor- 
relation function from ultrasoft gluon exchange, which appears first at NNNLO in the 
non-relativistic expansion, and is required for accurate top and bottom quark mass de- 
terminations from the threshold pair production cross section. The correction turns out 
to be large even for top quarks, but a definite conclusion on the attainable theoretical 
precision can only be drawn once the full NNNLO result, including potential and hard 
corrections, is available. A discussion of the sum of all known NNNLO terms for top 
quark production can be found in [TT] . 
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